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ABSTRACT 

Deep space observations of the JWST have revealed that the structure and masses of very early Universe galaxies at high redshifts 
(z ~ 15), existing at ~0.3 Gyr after the Big Bang, may be as evolved as the galaxies in existence for ~ 10 Gyr. The JWST 
findings are thus in strong tension with the ACDM cosmological model. While tired light (TL) models have been shown to 
comply with the JWST angular galaxy size data, they cannot satisfactorily explain isotropy of the cosmic microwave background 
(CMB) observations or fit the supernovae distance modulus versus redshift data well. We have developed hybrid models that 
include the tired light concept in the expanding universe. The hybrid ACDM model fits the supernovae type 1a data well but not 
the JWST observations. We present a model with covarying coupling constants (CCC), starting from the modified FLRW metric 
and resulting Einstein and Friedmann equations, and a CCC + TL hybrid model. They fit the Pantheon + data admirably, and 
the CCC + TL model is compliant with the JWST observations. It stretches the age of the Universe to 26.7 Gyr with 5.8 Gyr 
at z = 10 and 3.5 Gyr at z = 20, giving enough time to form massive galaxies. It thus resolves the ‘impossible early galaxy’ 
problem without requiring the existence of primordial black hole seeds or modified power spectrum, rapid formation of massive 
population III stars, and super Eddington accretion rates. One could infer the CCC model as an extension of the ACDM model 


with a dynamic cosmological constant. 


Key words: galaxies: high-redshift - cosmology: early Universe, dark energy, cosmological parameters. 


1. INTRODUCTION 


Observations with the JWST have revealed the existence of massive, 
bright galaxies in the very young Universe of age ^ 500 Myr (e.g. 
Naidu et al. 2022a, b; Curtis-Lake et al. 2023; Hainline et al. 2023; 
Labbe et al. 2023; Robertson et al. 2023), which is a small fraction 
of its current age estimate of 13.7 Gyr according to the standard 
ACDM model. Angular diameters of many such galaxies are an 
order of magnitude smaller than expected from the ACDM model 
(e.g. Adams et al. 2023; Atek et al. 2022; Chen et al. 2023; Donnan 
et al. 2023; Finkelstein et al. 2022; Naidu et al. 2022a, b; Ono et 
al. 2022; Tacchella et al. 2023; Wu et al. 2023; Yang et al. 2022; 
Austin et al. 2023; Baggen et al. 2023). Astronomers first identified 
such an ‘impossible early galaxy’ problem from observations with 
the Hubble Space Telescope (HST) at high redshifts z ~ 10 (Melia, 
2014, 2020, 2023). Dekel et al. (2023) stated ‘JWST observations 
reveal a surprising excess of luminous galaxies at z ~ 10’. According 
to Boyett et al. (2023) 'JWST observations confirm the existence of 
galaxies as early as 300 Myr and at a higher number density than 
expected based on galaxy formation models and HST observations'. 
Looser et al. (2023) observed the existence of a quiescent galaxy 
when the Universe was only 700 Myr old (see also Long et al. 2023). 
Bunker et al. (2023) wrote *Our NIRSpec spectroscopy confirms that 
GN-z11 is a remarkable galaxy with extreme properties seen 430 Myr 
after the Big Bang' (see also Tacchella et al. 2023); NIRSpec refers 
to ‘near infrared spectrograph’ and GN-z11 nomenclature means 
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Good North survey of galaxy with z = 11. Following spectroscopic 
confirmation of several photometric redshifts from JWST early 
galaxy observations, Haro et al. (2023) stated ‘our results solidifies 
photometric evidence for a high space density of bright galaxies at 
z > 8 compared to theoretical model predictions’. Related to the 
massive quasar analysis observed at z > 6, Eilers et al. (2023) wrote 
*.this quasar hosts a ten billion solar mass black hole less than 
1 Gyr after the Big Bang, which is challenging to explain with 
current black hole formation models’. ‘And the findings have been 
dazzling astronomers, revealing that stars and galaxies were forming 
and evolving much earlier than anyone had expected' as expressed 
by Alexandra Witze (2023) in a recent news article in Nature. How 
come the galaxies in the very early Universe were as evolved as 
those with billions of years of evolution, some as early as less 
than ~ 300 Myr after the big bang? These astonishing observations 
make the problem even more acute to resolve by tweaking well- 
established galaxy formation and cosmological models developed 
to satisfy lower redshifts observations (e.g. Haslbauer et al. 2022; 
Inayoshi et al. 2022; Kannan et al. 2022; Keller et al. 2022; Regan 
2023; Yajima et al. 2022; Atek et al. 2023; Mason et al. 2023; 
McCaffrey et al. 2023; Mirocha & Furlanetto 2023; Whitler et 
al. 2023a, b). Attempts have been made to compress time for the 
formation of population III stars and galaxies more and more, such 
as by considering the presence of primordial massive black hole 
seeds, and super-Eddington accretion rates in the early Universe 
(Ellis 2022; Bastian et al. 2023; Brummel-Smith 2023; Chantavat 
et al. 2023; Dolgov 2023; Larson et al. 2023; Maiolino et al. 
2023). While analysing GN-Z11 JWST-NIRSpec data, Maiolino et 
al. (2023) concluded that the black hole seed of this exceptionally 
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luminous galaxy at z = 10.6 must be accreting at an episodic rate of 
about five time the Eddington rate for 100 Myr (since z ~ 12 — 15) 
and is challenging for theoretical models (see also Schneider et al. 
2023). Chen et al. (2023) considered the presence of massive dark 
matter halos (see also Mauerhofer and Dayal 2023). As discussed 
elegantly by Melia (2023), both are considered unrealistic based 
on the zero angular momentum argument and observations showing 
most of the distant quasars accreting below or at the Eddington limit 
(Melia 2023). 

As stated by Wang and Liu (2023) ‘JWST high redshift galaxy 
observations have a strong tension with Planck CMB measurements'. 
They could not resolve the tension using alternative cosmological 
models either, including dark matter—baryon interaction, f(R) gravity, 
and dynamical dark energy (see also Santini et al. 2023). Parashari 
and Laha (2023) have shown that a blue-tilted power spectrum 
could potentially alleviate the tension with a low to moderate star 
formation efficiency. Studying the Balmer breaks at highest redshift 
galaxies, Steinhardt et al (2023) proposed it as a test of the ACDM 
model. They emphasized that the existence of stronger Balmer breaks 
out to z 2 11, will either demonstrate the early galaxy formation 
templates are invalid at high redshift or imply new physics beyond 
the ‘vanill? ACDM model. Lovyagin et al. (2022) have shown that 
the impossible early galaxy problem can be resolved amicably with 
the tired light theory in a steady-state Universe, first advocated 
by Zwicky (1929) to explain the early redshift observations by 
Hubble (1929). They have succinctly reviewed the current status 
of the problem and suggested some solutions. However, the tired 
light model cannot explain the extreme directional uniformity of the 
observed cosmic microwave background (CMB) radiation (Penzias 
and Wilson 1965). Additionally, the standard tired light model does 
not fit the supernovae type la data except at very low redshifts. 
An expanding Universe model can easily account for the observed 
redshift of distant galaxies and the CMB isotropy but has problems 
with the early Universe observations. However, the two approaches 
are not mutually exclusive. Gupta (2018a, b) has contemplated the 
existence of tired light in an expanding Universe. Our main objective 
in this paper is to explore if hybrid tired light and expanding 
Universe models can resolve the impossible early galaxy problem 
without conflicting with the relatively low redshift z < 2.5 data from 
supernovae type 1a (SNe Ia) observations in distant galaxies such as 
the Pantheon 4- data (Scolnic et al. 2022). Since the distance travelled 
by light reaching us is the same irrespective of whether in the tired 
light scenario or the expanding Universe, it constrains the parameters 
in them; one does not need an extra parameter to include tired light in 
an expanding Universe model. We will use two expanding Universe 
models to develop hybrid models, compare them, and see how they fit 
the SNe Ia and the galaxy size data from HST and JWST observations. 
The first hybrid model combines the ACDM model with the tired 
light model and is named the ACDM 4- TL model. The second hybrid 
model incorporates a new model derived from covarying coupling 
constant (CCC) approach (Gupta 2022). We start by defining an 
appropriate metric for the CCC approach and develop Einstein and 
Friedmann equations for the new CCC cosmological model; the 
hybrid model is dubbed the CCC + TL model. 

Section 2 of this paper develops the theory of the expanding 
Universe models and their hybrid counterparts, especially for the 
CCC model, and shows how it relates to the standard ACDM model. 
In Section 3, we attempt to fit the Pantheon 4- data (Brout et al. 
2022; Scolnic et al. 2022) for the models developed in Section 2 
to qualify them for testing further with the HST and JWST data. 
We discuss our findings in Section 4 and present the conclusions in 
Section 5. 
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2. THEORY 


2.1. Basics 


We will derive the Friedmann equations starting from the Einstein 
equations and see how they are modified when G and c are 
varying such that G ~ c? and the speed of light is used to measure 
distances (Gupta 2022). The basic Einstein equations (without the 
cosmological constant) are 

| 8zG (t) 

"C ey 

where G is the Einstein tensor and T,» is the stress-energy tensor. 
The FLRW metric takes the form 


ds? = d? f! — ary fy 


dr? 2 2 $29 2 
x (Sat (d0? + sin oq). (2) 


Tyv, (1) 


where f (t) = exp(a(t — to) defines the variation of c and the 
distances are measured using the speed of light; fo is the current 
time when f (t) = 1, a is the scale factor, and k (—1, 0, 1) 
determines the geometry of the space. The most significant difference 
from the standard FLRW metric is that the goo metric coefficient 
is time dependent (Gomide & Oehara 1981). Since f/f = o 


andf /f =a? , and G,, is a diagonal matrix, the essential 
components of the Einstein tensor are 
PEL , à fàM 
Goo — 3c — to t2a-t4|[- (3) 
a a a 
G i P TUN TAS NIS 
= a o — -= -]. 
" c? (1—kr?) (a a a a 


The stress-energy tensor for a homogeneous and isotropic universe 
considered a perfect fluid, has the form 


Thy = c? (€ + p) Uyty — pgu»- (5) 
Here & is the energy density, p is the pressure, g,,, are the metric 
elements of the FLRW metric, and the four-velocity u, are related 
through g""u,, u, = c?.Fora comoving observer, u,, = (c, 0, 0, 0). 


Since C, satisfies the contracted Bianchi identities and T,,, obeys 
the local conservation laws, we have 


82G (t) 
4 Tuv 
c(t) 
Since T,, is also a diagonal matrix, and spatial components are 


all equal, we need to consider only Too and Tj;. Now goo = f? , 
g% = f^? . Therefore, gu, uy = c? = uo uo = f?c? , and 


V^ Guy = 0, i.e. V" ( ) = 0, and V^ T,, = 0. (6) 


To =(e+ p P -pP =ef (7) 
_ pf^a 
T= I pZ (8) 


Now, £ and p both are dimensionally the same. Also, £ = pc? ~ 
mc? [r? ~ f?/ f? ~ f7! when m does not vary with time. It means 
Too and Tj, both scale as f.Sincec ^ f and G ~ f? 8x GT, Jc ~ 
$5. Therefore, V^ (81 GT, /c^) = 0 when V" T,, = 0 with the 
standard definition of 7,,, without the function f?. Einstein equations 
corresponding to the metric of equation (2) may now be written 


C a a IT 

kc 04 > 81G 

=z tea + 2a- + = me (9) 
a a a 3c 

kc? a åN? ä 8xG 

— o + 4a + 4+2-= ; P (10) 
a a a a C 
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Rearranging them leads to the Friedmann equations 


2NA p) $ 

å 8r G kc 2 a 
(=) CA ae (« x) bd 
: 4 
NRI NE EM (12) 
a 3c? a 


From these equations, we obtain the continuity equation 
a 

&-3— (e+ p)= —a (e 3p). (13) 
a 


The solution of this equation for the matter-dominant (p — 0) and 
radiation-dominant (p = 2/3) epochs of the Universe are, respec- 
tively, 


e = soa ^ exp (—a (t — to)) = eoa? fl, and 
e = soa exp(—2a (t — to) = £o af: (14) 


We will label these equations as representing the CCC (covarying 
coupling constants) universe. Comparing equations (11) to (14) with 
corresponding equations for the ACDM model, 


LAD 2 

a 8r G kc A 

-— 1 
(5) 3c 5 0g T 3' E 
Eu oe BD (16) 
a 3c? d 3" 
i43 (e+ p) 2 0, (17) 
a 

e = ea, ande = ena“, (18) 


immediately reveals: (i) the cosmological constant A of the ACDM 
model is replaced with the constant œ in the CCC model, (ii) the 
continuity equation has an additional term involving o in the CCC 
model, and (iii) energy density evolution has extra factors that must 
be considered for the CCC model. 


2.2. CCC model 


Defining the Hubble expansion parameter as H = à/a , we may 
write equation (11) for a flat universe (k — 0) as 


8rG 8rG 
Grale ied qma e (19) 
3c 3c 


In the matter-dominated Universe of interest to us in this work, 
using equation (14), we may write 


H = (Ho  o)a P f 0/9 — a. (20) 


Since the observations are made using redshift z, we have to see 
how the scale factor a relates to z in the CCC model. Along the spatial 
geodesic (0 and $ constant) between the observer and the source at 
a fixed time f, using the modified FLRW metric (equation 2) 


ds = a (t) f (t)dr. (21) 
Thus, the proper distance for co-moving coordinate r [recall that 
a (tg) = 1= f(to)l 

d, — a(t) f (t) fdr = a(t) f (r 2 dp (to) = r. (22) 


Since the light follows the null geodesic, equation (2) for a light 
emitted by a source at a time t, and detected by the observer at a time 
to yields 


to dt t 
bcc dard = d, (to). (23) 


On JWST early Universe observations 3387 


Thus, the expression for the proper distance is the same when using 
the standard FLRW metric. It can now be shown (Ryden 2017) that 
a = 1/(1 + z), ie. the same as for the ACDM model. 

The next thing to consider is to transpose f(t) to f(z), as it is the 
latter that we will require in calculating the proper distance. We may 
write equation (20) 


d 
t L eig + (Ho + a) f£ UBI = wa + (Eo + a)a O0 


dt 
x exp (-“S*) : (24) 


Its analytic solution (using WolframAlpha), with the boundary 
condition a = latt =f, is 


(Ga ( ae) ( a) 
a= exp +1 
2 aq 2 2 «a 


i 2/3 
x exp (=) : (25) 


2 


It can be written as a cubic equation 


3 (Hi 
Ax? - Cx - D— 0, withA— 1 E "Ru D C, 
a 
F=f 
D = —a?", andx = exp (-“ 2 2). (26) 


Its solution is 


f^-x- 


D DM cN 

Z+(( =) +(§)) 
D DW? c\3\'" 

my (( =) +(§)) (e 


Since the scale factor a — 1/(1 4- z), we have D — —[1/(1 4- ap^. 
Thus, the above equation provides the function f(z, Ho, o). 

Next, to determine the luminosity distance dz of an object, we 
need to understand how the photon energy flux evolves with the 
redshift in the CCC model. The photon energy is reduced by a 
factor (1+ z). The effect of the interval between photo arrival 
time on the flux is an additional factor (1 + z). Since in CCC, the 
distance scales as the speed of light, their evolution cancels out: two 
emitted photons separated in time by ôte are separated in space by a 
distance ôr, = c, ôte. This distance becomes ôro = Co ót,(1 + z) and 
the corresponding time interval to = codt.(1 + z)/co = ôte(1 + z). 
Thus the luminosity distance scaling is the same as for the ACDM 
model, i.e. D; = d, (to)(1 + z), and the expression for the distance 
modulus jz also the same, viz., 


u = 5log, (d; / 1Mpc) + Slog), (1 + z) + 25. (28) 


We will now focus on determining the expression for the proper 
distance d,. Since dt = dt x da/da = da/a, we may write equation 
(23) 


ides Sei c a Q9) 
p — Cc = E = _ ic , 
inia w a(t) a, ad ue a° H 
and since a = 1/(1 +z), da = — dz/ (a +z} = —dza?, we 


get, using equation (20) 
dz 
(Hy + a) 1 + 2) STP — a 


zd 
dy() = c f 5 — (30) 


own 
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2.3. ACDM model 


The corresponding expression for the ACDM model (matter domi- 
nant, flat universe) is (Ryden 2017), 

z dz 

S ; (31) 
9 4/Qmo( + 2° +1 = mo 


dp (to) = £ 

p V0) = Ho 
where Q2,,,0 is the relative matter energy density and (1 — €2,,9) is 
the relative dark energy density. Notice that both the models have 


only two unknowns, i.e. Hp and a for the CCC model, and Hp and 
Qm, for the ACDM model. 


2.4. Tired light model 


We would also like to consider the tired light model as it fits 
admirably well the high redshift galaxy size data (Lovyagin et al. 
2022), especially from the JWST. In this approach, one assumes that 
the photon energy loss, Adv, is proportional to its energy, hv, and the 
distance travelled, dr, i.e. dv = K vdr, with K as the proportionality 
constant. Thus, 


d vo d 0 
= Kare} =k far n (2) 
v be V d, Ve 


Ao 
= —Kd, > ln KJ” Kd,. (32) 
Here v, and v, are the emitted and observed photon frequencies and 
Ae and Ao are the respective wavelengths. The redshift is defined 
as z = (Ao — Ae)/ ^c, i.e. (1 + z) = A^o/ A;. Therefore, for the tired 
light 


1 
dp In +2). (33) 


It must reduce to Hubble’s law d = cz/Hp in the limit of z < 1, 
yielding K = Ho/c. Thus, the proper distance in the tired light (TL) 
model is given by 


C 
dy =F Ind 2). (34) 


2.5. Hybrid models 


Since the distance travelled is the same in any model, the parameters 
of a hybrid tired light and expanding universe model are obtained by 
equating d, for the two. Thus for the hybrid CCC + TL model, we 
equate equations (30) and (34), and write (with subscript t for tired 
light and c for the CCC model) 


Zc dz C 
c = In(1 + z;). 35 
S Het and +28? FQ —a mont s 
We will now determine H. and H, by considering the above 
expression in the limit of very low redshifts (z «& 1) and keeping 
the terms with up to second order in z. Since the left-hand side 
involves an integral, we need to retain only up to first-order terms in 
the integrand. 

Let us first consider f(z). It is not practical to determine it 
from equation (27) in the small redshift limit z = 0. We go back 
to equation (24) and rewrite it with a — 1/(1 4- z): 

d 

2 =a (1 +z) — (Ho + æ) f + 2). (36) 
Since f~“/”) varies slowly as compared to (1 + z)*/?, we can assume 
it is a constant equal to unity (i.e. its value at z — 0) in the above 
equation for solving it. With boundary condition t = f at z = 0, the 
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solution in the limit of z « 1 up to the first term in z can be easily 
determined: 


dz 
dt 


a (1-2- z) — (Ho 4- a) (re 32) 


1 
-Ho — 5 (SH + 30) z = A + Bz, (37) 


with solution 


B B/a 
exp [B (t — to)] (1 + az) = (exp [a (t — 19)]) 


B B a/B 
(red) = s= (12) (38) 


Since A = —H, for the CCC case, the above expression for f reduces 
to 

a 
f= (1-52). (39) 


The left-hand side of equation (35): 

: -(/2) -1 
T (3/2) az 

c f dz (H. 4- à) (1 +z) 1 — — —a 
0 


"n 


E 3 az E 
= ejas. +a) (1 + x) (1 + =) -« 


Ze 3 
= ejas + CH, + o) (5 + 
0 
C à (H. + a) a 
dz |1 
A | 2H, (3 T =) | 


CZc (H: + a) a 
Z (1 " (3+ x) æ): (40) 


Thus, by expanding the right-hand side, we may write equation (35) 
as 


CZc (H. + a) a CZi Zi 

1 3 Lec (1 = =) 4l 
H, ( 4H, ( *x)«) H, 2 (H) 
Equating the first term on both sides of this equation and compar- 
ing it with Hubble’s law yields z./ H. = z;/ H; = z/ Ho. Since the 


cumulative redshift z obeys (1 + z) = (1 + z,)(1 + z+), in the small 
redshift limit, we get z = Ze + z;. Therefore, 


£o a Hr > Hy = HET NP TE 
HA, H, Ho H, Ze 


H, 
H. (1 4 x) = H, eH. (42) 


Equating the second term on both sides of equation (41) 


(H: + a) Qa Zt H, 
3 = = Ze => H, 
AH, ( T x) d um cd 


"ORT (345), (43) 


a H. 


Let us now consider equation (35). With the above findings, we may 
write 
j i 
c 
0 (Ae + al + 20/? f(z) 0/2 — g 


_ [H +a) aN] T 34 
«sees (og s Ira pee a» 


Given H. and a, this equation determines (1 + ze) for any (1 + z) 
value and, therefore, also (1 + z,). We can now use either of equa- 
tion (31) or equation (34) to determine the proper distance dp 
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Figure 1. Pantheon+ (Brout et al. 2022; Scolnic et al. 2022) data fit for the four models discussed in the text.Okay 


while replacing Hy and z in there with H, and ze or H, and 
z, as applicable. When determining the luminosity distance, only 
(1 + ze) is involved in time dilation, not (1 + z+). The luminosity 
distance in the hybrid model is then d; = d,(1 + z.(1 + z,)'”. 
Since (1 + z;) = (1 + z)/(1 + ze), the distance modulus for a hybrid 
model comprising the tired light and the CCC model becomes 


u = 5log;o (d, / 1Mpo) + 2.5log,, [(1 + z) (1 + z-)] +25. (45) 


We can follow a similar approach for the ACDM + TL model 
(with subscript / for the ACDM model). Equating equation (31) 
and equation (34) 


CU dz 


[6 
J = T Ind + z). (46) 
A : V 9o + zy + 1— Qm,0 H, 


In the limit of z; «& 1 and z, « 1, and keeping the terms up to second 
order in redshift, we get the same results as for the CCC model by 
comparing the first-order terms, but H, = 1.5, 9 H; by comparing 
the second-order terms, and equation (44) is replaced with 


C i dz 


" 
Hi 9 TT + zy + 1— Q m,0 


Inl--z). (47) 


_ C 
^ L5Q,uo0Hi 


3. RESULTS 


3.1. Fitting pantheon + data 


We will now attempt fitting the Pantheon + data (Brout et al. 
2022; Scolnic et al. 2022) for supernovae type la (SNela). We have 
considered the following six models in matter dominated Euclidean 
universe: 


(i) The ACDM model. 

(ii) The CCC model. 

(iii) The hybrid model, comprising the ACDM model and the tired 
light model labelled ACDM 4- TL. 

(iv) The hybrid model, comprising the CCC model and the tired 
light model labelled CCC + TL. 


(v) The tired light model TL, modified to include an unknown flux 
loss correction term, labelled TL+, added to the luminosity distance, 
i.e. 5 log (1 + z)’, with £ a data fit parameter. The luminosity distance 
was corrected only for the photon energy loss since the time dilation 
is not relevant for it. 

(vi) The standard tired light model, TL, with luminosity distance 
again corrected only for the photon energy loss since the time dilation 
is not relevant for it. 


We used the standard Curve Fitting tool in Matlab for this 
purpose by minimizing x? (Gupta 2020). The first five have two 
free parameters; thus, the degrees of freedom are identical. They 
are all within the error bars of the data (Fig. 1; only the first four 
model curves are plotted). The x? values of the first four models 
are within 0.1 percent of each other at 745.7 as its mean value 
(Table 1). Even the TL + model’s x? value is within 0.6 percent 
of the mean value. However, as could be expected, the x? value 
for the standard TL model is significantly larger at 2580. Thus, 
the Pantheon + data rejects only the TL model but none of the 
others. 

It should be mentioned that we have not taken into consideration 
the systematic uncertainties associated with the data. However, they 
would affect the fitting of all the models and x? values similarly and 
thus are not expected to impact the comparative study of the models. 

The models’ fit parameters and x? values, along with the de- 
celeration parameter and the age of the Universe, are presented in 
Table 1. It is relevant at this point to refer to the work of Lopez- 
Corredoira and Calvo-Torel ( 2022) wherein they compare the fits of 
several alternative models without dark energy by fitting Pantheon 
data (Scolnic et al. 2018). 

Unlike in a previous paper (Gupta 2020), wherein absolute mag- 
nitude Mg = —19.35 was used to convert the apparent magnitudes 
reported in the data to the distance moduli u, Pantheon+ data 
(Scolnic 2022) includes distance modulus, apparently based on 
{M}-_B = - 19.253. Since Ho is degenerate with Mg, the choice of 
Mg affects the determination of Hp for different models. However, 
all the models yield about the same value of Ho. Thus, we do not 
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Table 1. Pantheon + data fit parameters, x? values, the deceleration parameter, and the age of the Universe. 


Parameter Unit LCDM CCC LCDM + TL CCC + TL TL+ TL 
Ho Kms^!Mpc'! 72.99+0.34 72.70 + 0.34 72.74 72.62 72.46 +0.34 64.42 + 0.36 
Hy Km s^! Mpc^! 72.99 72.7 60.48 + 1.06 59.51 + 1.06 0 0 
Hı Km s^! Mpc^! 0 0 12.26 13.11 72.46 64.42 
Qm.0 NA 0.3508 + 0.0243 NA 0.1351 + 0.0109 NA NA NA 
a /Hy NA NA —0.4953 + NA —0.7997 + NA NA 
0.0246 0.0143 
B NA NA NA NA NA 0.6418 + 0.0196 NA 
x NA 745.4 745.6 745.3 746.5 749.9 2580 
Fit Params. NA 2 2 2 2 2 1 
DOF NA 1699 1699 1699 1699 1699 1700 
qo NA —0.474 —0.371 -0.797 —0.780 NA NA 
to Gyr 13.75 13.69 19.25 26.70 NA NA 


expect the choice of the absolute magnitude to affect the findings of 
this paper. 


3.2. Angular-diameter angle and distance 


These are essential for testing a model based on the data on the size 
of galaxies at high redshift, especially those recently observed with 
the JWST and discussed extensively in the literature (e.g. Lovyagin 
et al. 2022 and references therein). 

The angular diameter distance d, is defined in terms of the physical 
size ôl of an object and its observed angular size 80 as d4 = ól/ 80. 
Using the metric (equation 2), the object at a location (r, @), i.e. 
dr = 0 and dó = 0, at time ft has a size given by 


ds? = a(t? f (t. r?d0? => ds = a (t) f (1) r80 = ôl. (48) 


Therefore, with r as the proper distance d,, the angular diameter 
distance becomes 


d, =a(t) f (d,. (49) 


We have to be cogent of the fact that the scale factor a in a hybrid 
model relates to the redshift z, due to expanding universe only, i.e. 
a = 1/(1 + zx), whereas an observer measures the total redshift z. 
The function f(t), of course, has no significance for the non-CCC 
models. In addition, a(t) is irrelevant for the tired light model. Thus, 
when we know d, for a model, we can immediately compute da, 
and, therefore, the size of an object from its observed angular size. 

Plots of angular-diameter distance for the five models are shown in 
Fig. 2. Two of the plots jump out: (i) the TL + plot has no maximum, 
but that can be expected as no expansion of the Universe is involved, 
and (ii) the CCC + TL plot has its peak value about four times 
higher than the ACDM value at z > 10. This dramatically affects 
the observed 50 and brightness of the distant galaxies. The plots are 
shown in Fig. 3 for 10 kpc objects against the backdrop of measured 
angular sizes of galaxies from multiple sources (the same as cited 
by Lovyagin et al. 2022; those points without errors in the original 
data are arbitrarily allowed a 5 percent error), including the latest 
JWST data. While all models are satisfactory at low redshifts, only 
TL + and CCC + TL are acceptable for high redshifts. 

It should be mentioned that the size evolution of galaxies, espe- 
cially at z — 10, is not well measured. In particular, the measured an- 
gular sizes are subject to the distribution of stellar populations within 
galaxies and their light or mass profiles (e.g. outskirts of galaxies 
are likely too faint to be detected). Also, we have not addressed the 
uncertainties related to baryonic processes and feedback mechanisms 
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that may potentially be effective in shaping galaxy sizes of observed 
galaxies. Additionally, not all redshifts determined initially via 
photometric colour selection, have been confirmed spectroscopically; 
some have been found spectroscopically at much lower redshift, e.g. 
a bright interloper at Zspec = 4.91 that was claimed as a photometric 
candidate at z ~ 16 (Harikane et al. 2023). 


3.2. Age of the Universe and redshift 


Galaxies at high redshifts, especially those observed by JWST, appear 
as evolved and as massive as those at lower redshifts. How can this 
be possible, considering that the age of the Universe at very high 
redshifts was less than half a billion years as per the ACDM model? 
This subject has been extensively discussed in the literature (e.g. 
Adams et al. 2023; Atek et al. 2022, Chen et al. 2023; Donnan et 
al. 2023; Finkelstein et al. 2022; Naidu et al. 2022a and 2022b; Ono 
et al. 2022; Tacchella et al. 2023; Wu et al. 2023; Yang et al. 2022; 
Labbe et al. 2023), so we will not delve on it.- Nevertheless, we will 
try to see how the Universe ages under different models. 

Let us first consider the CCC model. As per equation (24), we may 
write 


da 
dt 


—aa + (H, +a) f 0/gq70/ >t 


p da 
= 1 —aa + (H; + a) f-Q/2g-0/2 (50) 


Here a, = 1 for the age of the Universe today at ze = 0. To find the 
age for a specific ze, we replace with ac, 1/(1 + ze) and integrate. 
For non-hybrid models, z replaces ze. 

When working with a CCC + TL model, we have to find first the 
expanding universe component z, of z as described under equation 
(44), and then use equation (50) to determine the age corresponding 
to the observed z. Similarly, we determine the ages for the ACDM 
and ACDM 4- TL models. The Universe age has no meaning for the 
TL and TL 4- models. 

Fig. 4 shows how the age of the Universe decreases with increasing 
redshift for the four models. Fig. 5 shows the age increment with 
redshift for the two hybrid models compared to the currently expected 
age for the standard ACDM model. Both the hybrid models show 
significant increases, but the CCC 4- TL model provides a 10 to 20- 
fold increase at redshifts 10 to 20, giving enough time (5.8 Gyr at 
z = 10 and 3.5 Gyr at z = 20) for large, massive galaxies to form. In 
contrast, for the ACDM + TL model, corresponding age increases 
are 1.7 and 0.7 Gyr. 
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Figure 2. Angular diameter distance according to different models plotted against observed redshift. 
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Figure 3. The angular size of 10 kpc objects for the five models against the backdrop of measured angular sizes of galaxies from multiple sources, including 


the latest /WST data and some pre-/WST data labelled as NJWST (provided by Lovyagin et al. 2022). 
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Figure 4. Age as a function of the observed redshift for the four models. 


MNRAS 524, 3385-3395 (2023) 


€z0z Ainr zz uo ysanb Aq eve Lzz//sgee/e/vzs/elonie/seuuujfuoo: dno oruepeoe//:sdiu Wo pepeojuwoq 


3392  R.Gupta 
100 
—— CCCHTL/LCDM 
= ——LCDM-TL/LCDM 
(o 
t 10 
(D 
60 
< 
1 


10 100 1000 


Figure 5. Age advantage of the two hybrid models over the ACDM model. 


3.3. Tired light contribution to redshift 


We were curious to know how the ratio of the tired light and 
expanding universe redshift evolves with the total redshift for the 
hybrid models. It is plotted in Fig. 6. For the CCC + TL model, the 
tired light fraction starts at 22 per cent at z ~% 0 and declines slowly 
to 20 percent at z — 10, to 18 percent at z — 20, and to 3 percent 
at z = 1000. For the ACDM + TL model, the corresponding values 
are 20 percent, 11 percent, 7 percent, and 0.3 percent. Despite their 
relatively small contribution, the tired light modifies the observed 
redshift considerably since redshifts are multiplicative through the 
relation (1 + z) = (1 + z,X1 +z). 


3.4. Deceleration parameter 


We will now calculate the deceleration parameter for the models. It 
is defined by the relation qo = —ä/(a Hà), Ho being the expanding 
Universe part of the Hubble constant at the current time. Let us 
first consider the CCC model. From the first Friedmann equation 
(equation 11) for a flat Universe 


Hg = a €) — (o? + 2a Ho) 
4nG > 
= E eo = (Hè + (a? + 2 Ho)) /2. (51) 


From the second Friedmann equation (equation 12) in the matter- 
dominated universe, using the above equation 


a Hy = z (Hj + (a? +2aHo)) 


H un LE (52) 
a> aa ^U alm" 


Similarly, qo for the ACDM model may be expressed (Ryden 2017), 
since Q4 9 = 1 — Q2, 9, as 
1 


qo = 2 2m.0 imd Qa.o =-1 + 1.5€2, o. (53) 


In the above expressions, we must use the expanding Universe 
component of the Hubble constant when applying them to the hybrid 
models. The deceleration parameters and the Universe's ages for 
different models are shown in Table 1. 
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4. DISCUSSION 


The main objective of this paperis to explore if the hybrid models that 
include tired light cosmology can explain the deep space high redshift 
observation of JWST on the large-scale structures of the Universe, that 
is, the structure and evolution of galaxies. These observations show 
that the very early Universe galaxies were almost as bright, massive, 
and structurally as evolved as the galaxies in the late Universe but 
with rather small angular diameters. Under the standard ACDM 
model, the physical sizes of these galaxies turn out to be about 
10 percent of their expected sizes. It is because, according to this 
model, the angular diameter distance of objects reaches a maximum 
of about — 1.7 Mpc at z — 1.6 and decreases at higher redshift 
leading to the angular diameter 0 of a standard size object increasing 
with increasing z. Thus, if one observes a smaller 8, itis interpreted as 
a smaller size object. So, if a model has no angular diameter distance 
maximum, such as for the tired light models, or if the maximum is 
shifted to higher redshifts and thus has a higher value, then observed 
smaller 0 will translate into larger object sizes. Observed over density 
of galaxies in the region of z ~ 10 (e.g. Whitler et al. 2023) can also 
be explained accordingly. 

However, before any model can be considered seriously, it must 
pass the most basic test fitting the supernovae type 1a data compiled 
as Pantheon + by Scolnic et al. (2022). Additionally, it should, 
in principle, not conflict with the observed isotropy of the CMB 
radiation. The standard tired light model TL does not pass the basic 
test, as is evident from its x? value of 2580 against all others less 
than 750 in Table 1. The modified tired light model TL + with 
added unknown flux loss term 5 log (1 + z)? in the distance modulus 
equation gives a good fit to Pantheon + with a x? value of 749.9 
when 6 = 0.6418 + 0.0196. The inclusion of the additional term 
indicates that the TL model is missing something and that missing 
something is possibly the expanding universe component. However, 
neither of the tired light models can explain the extreme isotropy of 
the CMB radiation. Nevertheless, we have included the TL + model 
in some of the comparative figures for discussion purposes. All the 
remaining models provide similarly great fits to Pantheon + data 
(Fig. 1 and Table 1) and involve expansion of the universe and thus 
are compliant with CMB isotropy. 

One may be concerned about how the peak luminosities of 
supernovae type la in the Pantheon + data will be affected due 
to the variation of G and other constants. The peak luminosity can 
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Figure 6. Tired light redshift as a fraction of the expanding Universe redshift. 
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Figure 7. Object size under different models as a multiple of its size using 
the ACDM model. 


be related to the Chandrasekhar mass Men (Arnett 1982; Garcia- 
Berro et al 1999; Gaztañaga et al. 2001; Wright & Li 2018) with 
Mua, ~ (hc/ GP. Since h ~ f?, c~ f, and G ~ f?, we find that 
Men ~ f°. Therefore, the coupling constants’ variation does not 
affect the peak supernovae luminosity in the CCC model.! 

Coming back to angular diameter distance, its peak value of 7.45 
Gpc at z — 12.5 is the highest among the four models of interest, 
and its profile is closest to the TL 4- profile (Fig. 2). Thus, a small 
Observed angular size 0 of a galaxy does not translate into a smaller 
galaxy size. This is depicted in Fig. 3, which shows how one would 
observe 0 of a 10 kpc object with increasing redshift under different 
models; the curves are superimposed over the 0 of many galaxies 
from JWST observations and some pre-/WST observations. We notice 
that the CCC + TL model curve is very close to the TL + curve, 
whereas the ACDM + TL model curve is not. This may mean that 
the ACDM + TL model (also, the CCC and ACDM models) do not 
faithfully represent the early Universe evolution. 

An alternative way to analyse the angular size data of the galaxies 
would be to see what each model yields as their physical size. 
Since the physical size is directly proportional to angular diameter 
distance, the size increase in different models compared to the 
ACDM model can be easily determined, as shown in Fig. 7. We 


!In an earlier paper (Gupta 2022), the Planck constant scaling was incorrectly 
taken ash ~ f rather than h ~ f?. 


see an insignificant difference using different models at redshifts 
less than 1. The difference becomes perceptible at higher and higher 
redshift. The multiplier for the CCC+TL model becoming 5.6 at 
z=5,9.5 at z = 10, 12.8 atz = 15, and so on. Astronomers became 
concerned about the galaxy sizes using the ACDM model only at 
z > 5, especially at z > 10. Consequently, they suggested models 
and their improvements to create large structures of the Universe 
in a shorter and shorter time (e.g. Adams et al. 2023; Atek et al. 
2022; Chen et al. 2023; Finkelstein et al. 2022; Naidu et al. 2022; 
Ono et al. 2022; Tacchella et al. 2023; Wu et al. 2023; Yang et al. 
2022; Labbe et al. 2023; Chantavat et al. 2023; Yung et al. 2023). 
Instead of attempting to compress the cosmic timeline for creating 
well-evolved massive galaxies in a very young Universe to resolve 
the impossible early galaxy problem with the ACDM model (Melia 
2023) and cosmological simulations (McCaffrey et al. 2023), it might 
be prudent to consider alternative models that stretch the timeline. 
An alternative scenario in support of the standard model has been 
proposed by Prada et al. (2023) by contending that during the early 
epochs of the Universe the stellar mass-to-light ratio could not have 
reached the values reported by Labbe et al. (2023). However, one 
would need to explain the existence of very high z quasars, such 
as those considered by Latif et al. (2023), and remarkably high IR 
luminosity of massive galaxies at z ~ 8 discovered by JWST (Akins 
et al. 2023). 

We have shown in Fig. 4 how the cosmic timeline is stretched 
in the two hybrid models. We show the age of the Universe under 
the two hybrid models relative to the age under the ACDM model. 
While the stretch is modest (1.40 times for the ACDM + TL model 
versus 1.94 times for the CCC + TL model) at the current epoch, it 
becomes very significant at high redshift (2.97, 3.56, and 4.16 times 
at z = 5, 10, and 20, respectively, for the ACDM + TL model, and 
7.76, 12.3, and 19.3 times for the same redshifts, respectively, for the 
CCC + TL model). The cosmic time stretch is thus up to 5.8 Gyr at 
z = 10 and 3.5 Gyr at z = 20; there is no need to invent new physics 
for the rapid formation of galaxies. 

It would appear that except for the CCC + TL model, none of the 
models are able to reasonably resolve the impossible early galaxy 
problem; the TL models are unrealistic as they cannot reproduce 
the isotropy of CMB, and others do not provide realistic sizes of 
the early galaxies. Thus, the only realistic option is the CCC + TL 
hybrid model. 

The CCC model is derived from Dirac’s hypothesis of varying 
physical constants (Dirac 1937). His analysis was based on his 
large number theory and predicted the variation of G and the fine 
structure constant a. As discussed by Uzan (2011), if one constant 
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varies, then others also must vary, and their variations correlated. 
From the local energy conservation laws, Gupta (2022) found that 
the variation of c, G, the Planck constant h, and the Boltzmann 
constant k must follow G ^ c? ~ h!5 ~ k!? when the distance is 
measured with c. This leads to G/G = 3ċ/c = 1.5h/h = 1.5k/k, 
and therefore their variation can be represented by a single dimen- 
sionless function f(t), i.e. if one constant varies, then all of them 
do. The standard procedure is to consider a power-law function or 
a linear function. However, we found the most convenient form 
to use is f(t) = exp[a(t — fo)] with œ determined by fitting the 
Pantheon + data as a = —0.66H; (= c/c); a may be considered 
to represent the strength of the coupling constants’ variation. Since 
G/G = 3¢/c, we get G/G = —2.0Hp, about the same predicted 
by Dirac from his large number hypothesis. While there have been 
multiple attempts to constrain |G/G| by various methods to several 
orders of magnitude lower values (e.g. Teller 1948; Morrison 1973; 
Sisterna & Vucetich 1990; Degl’Innocenti et al. 1995; Thorsett 
1996; Corsico et al. 2013; Sahini & Shtanov 2014; Ooba et al. 
2017; Genova et al. 2018; Hofmann & Miiller 2018; Wright & Li 
2018; Bellinger & Christensen-Dalsgaard 2019; Zhu et al. 2019; 
Alvey et al. 2020; Vijaykumar et al. 2021), they all consider other 
constants to be pegged to their current value. However, keeping 
any of the constants fixed automatically forces o = 0, and therefore 
G/G — 0. 

We believe the two-parameter (Ho, œ) CCC model is possibly an 
extension of the ACDM model. Examining the Friedmann equations, 
one can imagine the CCC model as the ACDM model with a 
dynamic cosmological constant: (i) A — A(t) = —3(o? + 2ad/a), 
consisting of a static term and a dynamic term in the first Friedmann 
equation, and (i1) only the dynamic term in the second Friedmann 
equation. In the CCC model, recalling that o numerically turns out 
to be negative, it replaces A as the source of Universe expansion. 
One can even define energy density corresponding to o? + 2ad/a. 
However, it remains to be seen how well the CCC model can fit 
the Planck data and explain other astrophysical and cosmological 
observations. 

It should be mentioned that by treating the goo metric coefficient 
time dependent in the FLRW metric (equation 2), Gomide & Uehara 
(1981) were able to show that the local inertial effects are dependent 
on the overall structure of the Universe; cosmological models with 
positive curvature are Machian, whereas open ones are not. The 
CCC + TL model could, thus, be used to test the applicability of 
Mach’s effect in the Universe. 

It remains to be seen if the CCC 4- TL model can help address 
or explain other cosmology problems. Initially, we would want to 
determine if the new model is satisfactory for explaining CMB, 
BBN (big-bang nucleosynthesis), and BAO (baryonic acoustic oscil- 
lations) observations. 


5. CONCLUSION 


JWST is perhaps playing the same role as HST did in the 1990s — 
reinventing cosmology. HST put the ACDM model on the pedestal. 
JWST is challenging standard ACDM. In this paper, we have 
attempted to show that an extension of the ACDM model with 
deemed dynamical cosmological constant, when hybridized with 
the tired light concept and parameterized with Pantheon 4- data, 
provides a model, dubbed CCC + TL, that is compliant with the 
deep space observation of JWST on the angular sizes of high redshift 
galaxies. It stretches the cosmic time, especially at high redshifts, 
to allow the formation of large galaxies. It eliminates the need for 
stretching and tuning existing models to produce such structures 
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in the early Universe, thus amicably resolving the impossible early 
galaxy problem. 
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